Dipolar Bose-Einstein condensate of Stationary-Light Dark-state Polaritons 
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We put forward and discuss in detail a scheme to achieve Bose-Einstein condensation of stationary- 
light dark-state polaritons with dipolar interaction. To this end we have introduced a diamond-like 
coupling scheme in a vapor of Rydberg atoms under the frozen gas approximation. To determine the 
system's dynamics we employ normal modes and identify the dark-state polariton corresponding to 
one of the modes. We show that in contrast to atomic dipolar ultra-cold vapors dark-state polariton 
Bose-Einstein condensates proposed here can be stable for a negative dipolar interaction constant. 

PACS numbers: 42.50.Gy, 03.75.Hh, 32.80.Rm 



Long before the breakthrough of Bose-Einstein con- 
densation (BEC) in atomic systems quasiparticles have 
been considered as promising candidates for condensa- 
tion at high temperatures Their high critical tem- 
perature, which is one of their advantages compared to 
the atomic approach, is due to its inverse proportion- 
ality to the constituents' mass. Despite early attempts 
the condensation of quasiparticles in form of micro-cavity 
polaritons and thin-film magnon has only recently been 
achieved @, These systems show, however, two major 
downsides: (a) a relative short lifetime of the quasiparti- 
cles with respect to the thermalization time and (b) the 
two-dimensionality of both systems. In particular (b) 
prevents the formation of true long-range order of the 
one-particle density matrix and only quasicondensation 
can be achieved. 

The lifetime of quasiparticles can be extended if 
one uses, as shown in |4|], yet an other quasiparticle: 
the stationary-light dark-state polariton (SL-DSP) @ a 
quantum superposition of photons and collective Raman 
excitations. Thermalization can be achieved by intro- 
ducing a Kerr nonlinearity in the medium which effec- 
tively leads to a contact interaction of SL-DSP Here 
we consider a different mechanism of interaction based 
on an optically thick ensemble of Rydberg atoms. The 
Rydberg atoms induce a non-local long-range interaction 
between SL-DSP and hence introduce a novel interaction 
between polaritonic quasiparticles not known from earlier 
approaches. This can be used to simulate the dynam- 
ics of dipolar ultra-cold atomic systems. Moreover we 
will show that by varying the orientation of the dipoles 
one can control the interaction potential and realize new 
many body phenomena, e.g. very stable 3D BEC of DSP 
with long-range interaction or break the radial symmetry 
of dipole-dipole interaction. 

Let us consider a medium consisting of an ensemble 
of 4-level atoms with a diamond configuration resonantly 
interacting with four fields (Fig. [T]). The weak probe 
fields E± couple to the transitions \g) <-> |e±) and strong 



laser fields are present in only one leg of the diamond 
level configuration, e.g. E + = f2_ = 0, the scheme re- 
duces to a 3-level ladder-configuration which is known to 
support electromagnetically induced transparency (EIT) 
Q . This prominent effect has also been observed in Ryd- 
berg systems @, H| . Alongside transparency a probe field 
experiences at two-photon resonance a dramatic change 
of the refractive index properties of the medium leading 
for instance to a slow group velocity of a probe field pulse 
|. 

In the case that counter-propagating coupling fields 
in both legs are switched on and have comparable in- 
tensities the probe fields form a quasi-stationary pattern 
usually termed stationary light 

Q SEMI- We here 

consider the dynamics of stationary probe fields E± in 
an ensemble of Rydberg atoms as of Fig. [TJ The prop- 
agation of the probe fields is determined by the exter- 
nal control fields f2c which we assume to be sufficiently 
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FIG. 1: (Color online) The diamond linkage pattern (a) shown 
together with (b) the used field configuration for the counter- 
propagating control (fi±) and induced probe fields (E±). \g) 
denotes ground and \r) Rydberg state. \e±) are intermediate 
states. The one-photon detuning of the probe field transitions 
is denoted by A. The dipolar interaction of the Rydberg 
atoms induces an effective two- photon detuning A_r_r. It is 
assumed that ui c ^ to and thus there is no cross-couplings 
between the probe and coupling fields. Phase-matching is 
achieved by using non-parallel probe and coupling fields. 



control fields Q± drive the transitions 
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strong such that the medium shows no influence on them. 
For the later discussion we introduce slowly varying (nor- 
malized) probe field amplitudes £±(r, t) using E±(t, t) = 

^J1^£±{y, t) cxp [— i (cut =F kz)\ and control field ampli- 
tudes fi±(i) via Qc = fi±(i) exp [— i (u c t — k c r)]. 
In addition we introduce collective spin-flip operator 

= EjeV(r) lA t )jj( iy l/ 7V "( r ) where \p)j> k)j denote sin- 
gle atom eigenstates and the averaging is taken over a 
small volume V(r) around the position r with N(r) atoms 
inside |13|. 

The scheme shown in Fig. [T] as well as other linkage 
patterns [B[ suffer from a coupling of the probe £± (con- 
trol n±) field to the transition |e T ) <-> |r) (|<?) <-> |e±)) if 
the transition frequencies of the lower and upper legs of 
the diamond structures are comparable. In one dimen- 
sional systems discussed so far [HI, [H|, comparable 
transition frequencies have been employed to ensure EIT 
in hot atomic vapors with Doppler broadening Q. How- 
ever, in order to avoid this cross-coupling we assume that 
frequencies of probe and coupling fields are not equal, i.e. 
lu c =/= uj. This can be achieved by using the setup shown 
in Fig. 1 (b). 

We note that an efficient 4-wave mixing process can 
only be achieved if the setup used fulfills a phase- 
matching condition. If we denote the wave vector com- 
ponent of the control field parallel to the propagation di- 
rection of the probe field by kj! ± and the perpendicular 
component by k^ ± it is easy to show that the condition 
Ak = k + + k c — (k_ + k C) _|_) = can be satisfied if 

one chooses kj! ± = k T = ±fce 2 and hence k;f , = k;f_. 
Due to the non- vanishing perpendicular wave vector com- 
ponents of the control fields, we are restricted to the 
regime of cold Rydberg vapors. The choice above also 
eliminates the decay of stationary-light [l4[ . 

The light-matter interaction Hamiltonian including 
dipolc-dipolc interaction reads after rotating wave ap- 
proximation and after eliminating spatially fast changing 
contributions in the direction of propagation of the probe 
fields 
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where N is the number of atoms, Vt denotes the inter- 
action volume, and A is the one-photon detuning of £± 
and e (r) = £/p 2 (l — 3 cos 2 </j)|r| -3 denotes the frequency 
shift of the Rydberg state due the dipolar interaction of 
the atoms. We assume that the dipole moments p r of 



the Rydberg atoms are induced by a polarizing exter- 
nal field and thus are perfectly aligned 15|. The an- 
gle between the polarization direction of the dipoles and 
the difference vector r' — r is denoted by ip. It is well 
known that by using a rotating polarizing external field 
it is possible to tune the dipole-dipole interaction, i.e. 
to reduce its effective strength U and to even change its 
sign flri - Hit . This feature will be of importance later on. 
In eq. {l} the atom-field coupling constants are given 

by 9 = P-J 2he V t w ^ n P being the dipole moment of 
the \g) O |e±)-transitions. Finally we note that without 
Rydber-Rydbcrg interaction term the two photon detun- 
ing in each leg is zero. 

In the following we describe the derivation of an effec- 
tive equation of motion for the stationary-light dark-state 
polariton corresponding to the considered diamond-like 
coupling scheme. The propagation of the probe field is 
described by the Maxwell equation for the slowly varying 
amplitude in paraxial approximation 



d . c 2 



£± = igNa 
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The transversal dynamics is described by the last term on 
the left-hand-side, where = [d x ,d y , 0] is the transver- 
sal vector differential operator. The dynamics of the 
atomic system is governed by the Heisenbcrg-Langevin 
equations 
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where 7 M „ are transversal decay rates and are 6- 
correlated Langevin noise operators. 

We are interested here in the weak probe field limit, 
which allows analytical considerations. In this limit we 
assume that the Rabi frequency g of the quantum fields 
E± are, much smaller than £l± and the number of pho- 
tons in the input pulse is much less than the number of 
atoms. In this case we can treat the Heisenberg-Langevin 
equations pcrturbatively in the parameter g{£±) /£l±. 
We assume that all atoms are in the ground state in ze- 
roth order, i.e. a gg = 1. In first order we find 

d 

-QjVge± (r, t) = -r se± <7 se± + ig£± + i^e lk - r a gr + F ge± 
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—a gr (r,t) = -iA RR (r,t)fr gr + ie~ lK r ^ &ge-jtoj 



(5) 



withA fliJ (r) = (N/V t ) J v rr (r', i)e(r'-r)dV. Note that 
the Rydberg state \r) is long-lived and that the decoher- 
ence rate j gr is negligibly small Q. It would lead to a 
small loss rate of the considered SL-BEC which we are 
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going to disregard here. Because of the above assump- 
tion a Langcvin force term F gr is not necessary. For 
the sake of simplicity we assume that 7 ffe± = 7, hence 
r ffe± = r = 7 + iA. Within the scope of the publica- 
tion we restrict ourselves to the case O = f2 + = il_ and 
also assume that we can choose the control field Rabi 
frequency real [l2| . This is the limit of pure stationary- 
light. Equations flD, © together with eq. ([2]) form a 
self-consistent set of equations and determine the dynam- 
ics of the system. We note that eqs. (j4]) and ((U can be 
recast in the form 



^5(r, t) = -TS + ig£ s + iy/Xl(t) 



j-Mr.t) = -rD + igE D +F D , 
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-<T sr (r, t) = -iA RR (r, t)a gr + iV2n(t)e~ ik c- S, (8) 

where we have introduced sum respectively difference po- 
larization via s/2S = <Jg e+ +o ge _ and y/2D = cr ge+ —o ge _ . 
In addition we have defined the new field modes via 
y/2£ s = £++£- and y/2£ D =£+-£- [H Within 
this new basis the ongoing dynamics becomes more trans- 
parent. 

For a sufficiently long probe field pulse with character- 
istic pulse length T we can adiabatically eliminate the 
difference polarization mode D using cq. ([2j and eq. 
if irrl > 1 and find \u\ 
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To obtain eq. (j9]) we have furthermore employed that the 
characteristic length of the probe pulse is much larger 
then y/\ iA/j + 1| L^/k where L abs = jc/g 2 N is the 
absorption length of the medium in absence of EIT. Using 
the above result (|9|) and the shortened wave equation for 
the sum mode ([2]) we find that the sum polarization is 
given through the sum mode by 
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(10) 

Next we consider eq. ([5]) in the limit of an adiabatically 
slow changing control field. This limit can be character- 
ized by two features: (a) the Langevin noise operators 
can be neglected [B| and (b) the coherence between the 
ground and the Rydbcrg level is given by 
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ij,,(Anii(r,t)) can easily be satisfied. Here (Aim(r,t)) 
is the averaged detuning caused by the dipolc-dipole in- 
teraction and Lp U i se is the length of the probe pulse in 
the medium [2^. Combining eqs. ((HJ, (fl0|) and the result 
for the spin coherence ([lip we arrive at 

d ( A\ d 2 c -\ fi , , 



g 2 N ( d 
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Finally we introduce the dark-state polariton operator of 
the diamond-type coupling scheme. We can identify the 
system's unique dark-state polariton using the procedure 
introduced in [j| 
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#(r,t) = cos9£ s (r,t) - sm6V Na gr {r,t)e 



where the mixing angle 6 is defined via tan 2 6 = 
g 2 N/2Q 2 . In the adiabatic limit equation eq. (fT3|) can be 
simplified to £s(r,i) = cos 6'5'(r, t). By substituting this 
expression in eq. (|12[) we finally arrive at a Schrddinger- 
like propagation equation for the dark-state polariton op- 
erator 
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where the dipolar frequency shift in the polariton basis 
has the form 6 RR = (N/V t ) sin 2 (<9) / *t( r ')$( r ') e ( r ' _ 
r)d 3 r'. In eq. (HH we have introduced the transversal 



(9) and longitudinal mass using m± = hk/v g 
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where a = (2kL abs (A/7 -i)) . The group velocity of 
the probe field in the EIT medium, which can be con- 
trolled by changing the intensity of the coupling field, is 
given by v gr = ccos 2 9. 

As noted by Fleischhauer and colleagues Bose-Einstein 
condensation of SL-DSP can be achieved by raising the 
critical temperature of the polariton gas using the vari- 
able effective mass of the quasiparticles Q ■ In contrast 
to [H our scheme allows for direct thermalization via the 
dipolar interaction of the polaritons. Well below the crit- 
ical temperature the system can be described using the 
order parameter 4>(r,t) = y/N/V t (9(r,t)) [2l|. Due to 
this eq. (|T4|) simplifies to 
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If we denote the characteristic time for changes of the 
control field Rabi frequency by T and if we assume 
that the density of excitations in the system is small 
then the adiabaticity conditions O 2 j '7 'a/ 'L a bs I 'L pu ise S> tuning is large compared to the relaxation rate (7 <C A) 



which describes the dynamics of the so called dipolar 
Bose-Einstein condensates without contact interaction 
12, 22, HI . Note, that the longitudinal and transversal 
masses in eq. (fi"5|) arc different. In general the longitu- 
dinal mass mn is imaginary, but, if the single photon de- 
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the imaginary part of mil can be neglected. Equation 
(fl"5|) is the main result of this manuscript and shows that 
the dynamics of purely dipolar BECs can be mimicked 
using stationary light in EIT media. 

One of the major problems of dipolar BEC are the 
instabilities caused by attractive forces between aligned 
dipoles. In order to study the stability properties of (fT5|) 
we calculate the dispersion relation of weak density per- 
turbations with frequency v and wave vector q of a ho- 
mogeneous dipolar SL-BEC with density Udsp (HJ and 
find 



2m 
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C dd (3cos 2 /3- 1) 
(16) 



where C d d = 8irUp 2 n d s P - Moreover, qj_ = ^jq 2 + q 2 re- 
spectively q z are the transversal and longitudinal com- 
ponents of q. The angle between the excitation wave 
vector q and the direction of dipoles we denote by /3. In 
contrast to the usual dipolar BEC there is no contact 
interaction term in the present system [l?! 23 1. In the 
following we discuss two cases, defined by the orientation 
of the considered dipoles with respect to the direction 
of propagation of the probe fields, to show the unique 
properties of a stationary-light dipolar BEC. 

(i) Longitudinal orientation of dipoles 
Let us first consider the case when the dipoles are ori- 
ented along the z-axis. Then cos 2 j3 = q 2 / (q\ + q 2 ) and 
the Bogoliubov spectrum (|16|) can be written as 
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(17) 

where q = {q x , q y , q z a}. Since a < 1 (the typical value 
of a is on the order of 10 -4 ) the cloud is very stable (un- 
stable) if Cdd < {Cdd > 0), the only unstable (stable) 
region corresponds to the waves propagating perpendic- 
ular (parallel) to z direction due to the small value of a. 
Note that the interaction parameter Cdd can have both 
positive as well as negative values [3, Gj| . 

(ii) Transversal orientation of dipoles 
For the second case we assume that the dipoles are ori- 
ented along the y-axis. Then the angle between dipoles 
and wave vector q is defined by cos 2 /3 = q 2 / (q\ + q 2 ) 
and for the excitation spectrum we find 



q 2 
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(18) 

Equation (|18[) indicates that the radial symmetry of in- 
teraction is broken and that the dispersion relation be- 
comes strongly anisotropic in 3D. Namely, for the case of 
small a and Cdd > the system is stable (unstable) for 
the excitations propagating parallel to y-axis (x-axis). 



In summary, we have discussed stationary-light in an 
ensemble of cold Rydbcrg atoms. To this end we intro- 
duced a diamond-like light-matter coupling scheme where 
the upper state is given by a Rydberg state. We could 
show that the unique polaritonic quasipartilce behaves 
in the adiabatic limit like a Schrodingcr particle with a 
purely dipolar inter-particle interaction. Moreover, we 
could show, by analyzing the Bogoliubov spectrum of a 
homogeneous dipolar BEC, that for a special choice of 
the dipolar interaction parameter the considered dipolar 
SL-BEC is, in contrast to usual dipolar BECs, very sta- 

. blc. Depending on the orientation of the Rydbcrg dipoles 
used, the SL-BEC shows also a very anisotropic behav- 

'ior. The presented scheme opens up new roads between 
quantum optics and condensed matter physics and puts 
polaritonic quasiparticles into new light. 
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